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We report low-temperature specific heat measurements for high-quahty single crystalline KFe2As2 
(Tc ~ 3.5 K). The investigated zero-field specific heat data yields an unusually large nominal Som- 
merfeld coefficient 7n = 94(3) mj/mol K'^ which is however affected by extrinsic contributions as 
evidenced by a sizable residual linear specific heat and various theoretical considerations including 
also an analysis of Kadowaki- Woods relations. Then KFe2As2 should be classified as a weak to in- 
termediately strong coupling superconductor with a total electron-boson coupling constant Atot ~ 1 
(including a calculated weak electron-phonon coupling constant of Aei-ph = 0.17). From specific 
heat and ac susceptibility studies in external magnetic fields the magnetic phase diagram has been 
constructed. We confirm the high anisotropy of the upper critical fields ^oIIc2{T) ranging from a 
factor of 5 near Tc to a slightly reduced value approaching T = for fields _B || ab and || c and 
show that their ratio F slightly exceeds the mass anisotropy of 4.35 derived from our fuU-relativistic 
LDA-band structure calculations. Its slight reduction when approaching T = is not a consequence 
of Pauli-limiting as in less perfect samples but point likely to a multiband effect. We also report 
irreversibility field data obtained from ac susceptibility measurements. The double-maximum in the 
T-dependence of its imaginary part for fields B || c indicates a peak-effect in the T-dependence of 
critical currents. 

PACS numbers: 74.25.-q, 74.25Bt, 74.25Ha 



I. INTRODUCTION 

Since the discovery of superconductivity (SC) in an 
electron doped LaFeAsO (La-1111) compound with a su- 
perconducting transition temperature Tc 26 K [l[ , iron 
pnictides arc of great interest in fundamental condensed 
matter physics due to their large variety of structural, 
magnetic and electronic properties. In order to under- 
stand the nature of superconductivity in Fe pnictides a 
huge amount of theoretical and experimental studies has 
been performed but nevertheless many questions have 
not been answered yet, such as the symmetry of the 
order parameter and the pairing mechanism as well as 
their relation to the magnetic properties. In this sit- 
uation, low-T specific heat and magnetic susceptibility 
measurements might be helpful since they provide in- 
sight into many-body physics via the renormalization of 
such physical quantities as: the Sommerfeld coefficient 
7ci (a measure of the renormalized density of states), 
the irreversibility field Hirr, the upper critical field Hc2^ 
its anisotropy etc. all arc important factors which affect 
the superconducting and the normal state properties as 
well. In particular, they can shed light on the Fermi 
surface topology and other relevant aspects of the elec- 
tronic structure. To address the role of magnetism in 
the formation of the superconducting state the heavily 
hole doped KFe2As2 is worth to be studied due to its 
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distinctive characteristics with respect to other stoichio- 
metric 122 and 1111 Fe-pnictide compounds: (i) there is 
no static magnetic ordering in the sense of an ordinary 
spin-density wave (SDW) or an orthorhombic structural 
transition (ii) Superconductivity occurs in rela- 

tively dirty samples near 2.8 K Q and increases in cleaner 
high-quality single crystals up 3.5 to 3.7 K (iii) Re- 
markably, no nesting of the Fermi surface was detected 
in contrast to e.g. Bao.6Ko.4Fe2As2 [1, Q. However, a 
neutron scattering study of heavily hole doped supercon- 
ducting KFe2 As2 revealed well-defined low-energy incom- 
mensurate spin fluctuations at [7r(l±2(5), 0] with 5 — 0.16 
01 . This is different from the previously observed com- 
mensurate antiferromagnetism (AFM) of electron-doped 
AFC2AS2 (A = Ba, Ca, or Sr) at low energies. Addition- 
ally, de Haas- van Alphen [1] and cyclotron resonance 
studies of KFe2As2 revealed a strong mass enhancement 
of the quasi-particles. Furthermore, it exhibits a very 
large anisotropy as compared with less hole-doped mem- 
bers of the 122-family and other Fe-pnictide and chal- 
gogenide supperconductors [l3l After naturally electron- 
ically more anisotropic 1111 and Tli_j,Rbj,Fei„i5Se2 with 
y ^ 0.4,(5 ~ 0.3 superconductors showing only slightly 
larger or comparable slope anisotropics of ^ 5 to 6 K- 
122 belongs to the most anisotropic pnictides [llj. A 
complete understanding of their critical field slopes near 
Tc is still missing due to the complex interplay of pair- 
breaking impurities and the symmetry of the supercon- 
ducting order parameter [T^] ■ 

The magnetic phase diagram of KFe2As2 has been 
studied via resistivity measurements on single crystals 
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[J], however, its determination using thermodynamic 
bulk techniques on single crystalline material is lack- 
ing up to now [13, [3l ■ In the context of recently dis- 
covered very high surface upper critical field with 
Hcz{T)/Hc2{T) - 4.4 for the external field || to the ab- 
plane in Ko.rsFei. 63862 [lEl is noteworthy, because for 
that direction the nucleation starts already much higher 
at iJc3 and resistivity and or/ac susceptibility measure- 
ments might in principle lead to a confusion between Hc2 
and 7Jc3 suggesting this way an overestimated anisotropy 
of the upper critical fields. Hence, for this geometry, 
specific heat studies of high-quality single crystals are 
mandatory to address their bulk anisotropy. 

Here we present low-T specific heat and ac magneti- 
zation studies on high-quality superconducting KFe2As2 
single crystals with a larger and a much higher resid- 
ual resistivity ratio as compared to the first single crys- 
tals used for an upper critical field study for KFe2As2 
by Terashima et al. Q, where a large anisotropy ratio 
of the upper critical fields as well as for the electric re- 
sistivity perpendicular and parallel to the a6-plane have 
been reported. The obtained data are analyzed within 
the framework of various theoretical approaches. In par- 
ticular, we found good agreememt with the electron mass 
anisotropy derived from density funtional theory (DFT) 
considering in- an out-of plane plasma frequencies and 
Fermi velocities. 



II. EXPERIMENTAL 

Single crystals of KFe2As2 have been grown using a 
self-flux method with K:Fe:As in the molar ratio of 1:4:4. 
All preparation steps like weighing, mixing, grinding and 
storage were carried out in an Ar-filled glove-box. As a 
first step the appropriate amounts of the precursor mate- 
rials FeAs and Fe2As were thoroughly grinded in an agate 
mortar, secondly the exact amount of weighed metallic 
K was deposited at the bottom of an alumina crucible, 
where on top of it the well grinded mixture is placed 
carefully, finally sealed in a niobium crucible. The sealed 
crucible assembly is placed in a vertical furnace, heated 
up to 1373 K and cooled down to 1023 K with a rate of 
2 K/hour. Finally the furnace is cooled very fast from 
1023 K to room temperature. All crystals are grown with 
layer-like morphology and they are quite easy to cleave 
along the ah plane. The quality of the grown single crys- 
tals was assessed by complementary techniques. Several 
samples were examined with a Scanning Electron Micro- 
scope (SEM Philips XL 30) equipped with an electron 
microprobe analyzer for a semi-quantitative elemental 
analysis using the energy dispersive x-ray (EDX) mode. 
The composition was estimated by averaging over several 
different points of the platelet-like single crystals and is 
found to be consistent and homogeneous with a 122 struc- 
ture within the instrumental error bars. Typical crystal 
sizes with a rectangular shape were about 1.2x0.5 mm^ 
and a thickness of 50 /xm along the c-axis. All crystals 



exhibit very similar Tc- values (see the insets of Figs. 2 
and 8) within the experimental error of about 0.25 K. 

Low-temperature specific heat and ac magnetization 
have been determined using a Physical Property Mea- 
surement System (PPMS from Quantum Design). The 
specific heat data were measured using a relaxation tech- 
nique. For the measurements at H \\ ah a, small copper 
block has been used to mount the sample on the specific 
heat puck. The heat capacity of the copper block was 
determined in a separate measurement and its value has 
been subtracted from the raw data of KFe2As2. 



III. RESULTS AND DISCUSSIONS 
A. ac magnetization measurements 

Figure [T] depicts the T-dependence of the volume ac 
susceptibilities (x't, and x" v) of our KFe2As2 single crys- 
tal. The measurements were done in an ac field with 
an amplitude Bac ~ 5 Oe, a frequency f = 1 kHz and 
DC fields up to 5 T parallel to the a&-plane as shown 
in Fig. [1^ and parallel to the c axis in Fig. [T}). Special 
care has been taken to correct the magnetization data for 
demagnetization effects, where the demagnetization fac- 
tor has been estimated based on the crystal dimensions 
[l6j . The ac susceptibility measurements can be used 
for an investigation of the fiux dynamics in superconduc- 
tors The imaginary part x"v is related with the 
energy dissipation in a sample and the real part x'y is 
related with the amount of screening. Both these func- 
tions depend on the ratio between skin depth Ss and the 
sample dimension L in the direction of the fiux penetra- 
tion. In the normal state Ss ~ {Pn/ f)^'^, where pn is the 
normal state resistivity and / is the frequency. In the su- 
perconducting state the skin depth Ss oc Xl if an external 
magnetic field is below the first critical field Hd, where 
Xl is the London penetration depth. For magnetic fields 
above Hd, Ss oc Lb, where Lb ^ Bad Jc is the Bean's 
penetration depth and Jc is the critical current density. 
In general, if L <C (5^ an ac field completely penetrates 
the sample and thus the susceptibility is small. In the 
case, if i » Js, most of the sample volume is screened, 
therefore, A-kx' v = — 1 and x" v — >■ 0. In accordance with 
this the ac susceptibility data measured at low T confirm 
the bulk superconductivity of our KFe2As2 single crys- 
tal (Fig. [1]). Tc ^ 3.6(1)K has been extracted from the 
bifurcation point between x' v and x" v susceptibilities as 
shown in the inset of Fig. [T}d. This point is related with a 
change in the linear resistivity due to the superconduct- 
ing transition. It can be also used for the determination 
of the T-dependence of the upper critical field Hc2 from 
ac susceptibility data measured at various DC fields [17 1 . 

At T < Tc the function x" v increases with decreasing 
temperature and some value of T = Tirr where (L ~ ^s) 
x" V has a maximum (see Fig. [1]). The above discussed 
Bean's approximation predicts that the T-dependence of 
the peak in x" v follows the T-dependence of the critical 
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FIG. 1: (Color online) (a) The T-dependence of the complex 
ac susceptibility components A-kx v and 47rx"ii of KFe2As2 
has been measured in an ac field with an amplitude 5 Oe and 
a frequency of 1 kHz upon warming in different DC magnetic 
fields after cooling in zero magnetic field with (a) B \\ ab 
and (b) B || c. The sharp superconducting transition with 
~ 100 % superconducting volume fraction indicates the bulk 
nature of superconductivity and the high quality of our crys- 
tal. The inset shows the criteria used to obtain Tc and Tirr', 
for details see text. 



current density Jc- Thus, one might relate Tirr with an 
irreversibility temperature and use its DC field depen- 
dence to obtain the irreversibility field Hirr [3. How- 
ever, we remind the reader that Hi^r defined tliis way is 
not the 'true' irreversibility field since by definition Hi^r 
is the field at which Jc = 0. From this point of view 
it is better to use dc magnetization data to obtain the 
irreversibility line. In general, we observed a rough agree- 
ment between Hirr obtained from dc and ac susceptibil- 
ity data but the large value and the strong T-dcpendcncc 
of the normal state dc susceptibility lead to a large un- 
certainty in determination of the Hi^r- Therefore, to ob- 
tain the irreversibility line we used ac susceptibility data. 
(The DC susceptibility data will be presented elsewhere 
[20|). Thus, with some cautions we relate the maximum 
in the T-dependcnce of x" v to Ti^r- It can be seen in 
Fig. [TJd that at non-zero -Bdc||c the single maximum in 



x"v splits into two features at Tli 
fore, for SdcHc we defined two different 'irreversibility' 
fields Hlirr and H2irr- The T-dependence of these fields 
is plotted in Fig. [S] 



B. Specific heat studies 

Figure 2 shows the T-dependence of the zero-field spe- 
cific heat measured down to 0.4 K. A clear sharp anomaly 
was observed near 3.5 K in agreement with the magneti- 
zation data. In order to determine the zero-field normal 
state Sommerfeld coefficient 7n, the specific heat can be 
plotted for T > Tc as Cp/T versus T^ following 



7„T + /33T3 + /^sT^ + 



(1) 



with and /Ja, as the nominal electronic and lat- 
tice coefficients, respectively. The obtained values for 
our KFe2As2 sample are 7„ = 94(3) mJ/mol and 
h = 0.79 mJ/mol-K^ and 13^ = 6.09 • 10~'*mJ/mol K^. 
Our 7n value compares very well with 7„ = 93 mJ/mol 
reported in Ref. [2l| . From the relation for the Debye 
temperature, Od = (12tt^RZ /SP^) ^Z^, where R is the 
molar gas constant and Z = 5 is the number of atoms 
per formula unit, wc obtain 9]j — 214 K. 

Notice the enhanced value of 7n KFe2As2 as compared 
to other stoichiometric and nonstoichiometric 122 com- 
pounds or to any other superconducting iron pnictide or 
chalcogenide to the best of our knowledge. The low val- 
ues of 7n for BaFe2As2 and SrFe2As2 are not surprising 
since large parts of its Fermi surface are gapped due to 
the well-known magnetic spin density wave (SDW) tran- 
sition at high temperatures. Hence, in these cases a com- 
parison with a hole doped system where the SDW tran- 
sitions are suppressed is more meaningful. For instance, 
for the closely related, nearly optimal hole-doped systems 
Bao.68Ko.32Fe2As2 (Tc = 38.5 K) [H, Bao.6Ko.4Fe2As2 
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FIG. 2: (Color online) The T-dependence of the zero-field 
specific heat of KFe2 As2 for 400 mK < T < 200 K. The inset 
shows Cp/T versus T of our OT and 9T data together with 
the zero-field of another sample (S2) shows the same Tc{H) 
behavior down to 1.8K. 
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S= 36.5 K) g or Bao.esKo.asFesAsz (T, = 29.4 K) 
, the Sommerfeld parameters 7„ = 50.0, 63.3 and 
57.5 mJ mol^^K~^, respectively, have been reported. In 
view of their much higher Tc-vahies, ascribed by the au- 
thors to an essential contribution from strong coupling 
corrections with A ~ 2 [22| for the electron-boson cou- 
pling constant (spin fluctuation mediated interband cou- 
pling) and a comparable bare value of 7b ^ 10 mol~^K~^ 
(see below) according to DFT-band structure calcula- 
tions the unusual large value for KFe2As2 reported above 
provides a surprising puzzle. However, it can be some- 
what reduced, if there is an essential extrinsic contri- 
bution to the system of intinerant charge charriers, e.g. 
due to defect states with low-energy excitations [lO] ■ To 
be consistent with that analysis we are forced to as- 
sume, that KFC2AS2 is not in the strong coupling limit 
which seems to be natural in view of its low Tc-value (to 
be discussed within the framework of Eliashberg-theory 
elsewhere). The general situation, independent on the 
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FIG. 3: (Color online) The relative extrinsic linear specific 
heat coefficient of KFe2As2 using the experimentally observed 
nominal 7n — 94 mJ/mol -K^ vs. the total electron-boson 
coupling constant Atot given by Eqs. (4,7) for various bare 
7b-values obtained from density of states as calculated by 
various DFT codes [22,25,26] (see text) and using a typical 
high-energy renormalization factor of 77 = 2.7 (see Eq. (4)). 
The data points show the results of simulations within single 
band d-wave Eliashberg theory to reproduce Tc — 3.5 K and 
a spectral density for spin fluctuations adopted from recent 
INS data [7] and including also a weak electron-phonon inter- 
action and a weak Coulomb pseudopotential fi* (details will 
be discussed elsewhere). 
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FIG. 4: (Color online) The normalized specific heat jump 
at Tc vs. the phenomenological total electron-boson coupling 
constant A and various DFT-derived bare linear specific heat 
values [22,25,26] Red and blue symbols: results using the cor- 
responding 7ci values from Fig. 3 extended by two strong 
coupling points for 7ci — 7n and the experimental value of 
Acei/Tc = 45.6 mJ/ mole -K^) 



strength of the electron-boson coupling regime and the 
symmetry of the order parameter, is plotted in Figs. 3 
and 4 using the calculated DFT- values 76 [H, [2^, [Ig as 
convenient bare values. Thereby the high-energy renor- 
malization 77 = 2.7 as derived from the calculated DFT 
plasma frequencies of 2.58 eV [H and 2.56 eV [M [13 
and and the expected experimental unscreened plasma 
frequency of about 1.55 eV have been taken into account 
[2^. There are at least two experimantal hints which 
clearly point to the existence of such an extrinsic subsys- 
tem which manifest itself in a substantial residual linear 
specific heat visible at very low-T at ambient fields (see 
Fig. 8) and in high fields of about 9 T where the super- 
conductivity is well suppressed (see inset of Fig. 2). From 
the latter one estimates 701 < 70 mJ/mol K^. In the next 
section we will provide theoretical arguments in favor of 
a significantly reduced intrinsic Sommerfeld term 701. 



5 



C. Theoretical estimates of the thermal mass 
enhancement 

1. Kadowaki- Woods analysis 

The weak coupling result can be understood at a qual- 
itative level also by analyzing the so called Kadowaki- 
Woods relation (KWR) kkwr = Ap/j^, where Ap de- 
scribes the T^-contribution to the resistivity at very low 
T: p{T) = po + ApT^ observed so far only in few very 
clean samples [l^ with extremely large residual resitiv- 
ity ratio p{300K) / p{5K) « 500 and 7v is the volumetric 
Sommerfeld coefScient. The latter is related to the usu- 
ally used molar quantity 70 in the present case with two 
KFC2AS2 units per unit cell by the relation 

7„ = 27o/A^aK 

where A'a denotes the Avogadro's number and Vu = 
204.2 A3 is the unit cell volume of KFe2As2. Then fol- 
lowing Hussey [2^, one has for the case of a (quasi)-2D 
system with cylindrical Fermi surface sheets 
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Notice the cancellation of many-body renormalizations 
on the r.h.s. and the different exponent -1.5 for the den- 
sity n as compared with oc n"^ within a similar expres- 
sion proposed recently and given here for comparison also 
for the 2D-case, only: [30l |: 
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Since in the stoichiometric case of K- 122 there is exactly 
one hole in the three bands (i.e. per formulae unit) which 
cross the Fermi energy one obtains for the corresponding 
electron density ?ici = 4.94 x 10^ 



V2nnc. = 2.07- 10 



-10^ 



for the 2D-cffectivc Fermi wave vector of electrons. In- 
serting our value of fcp.oi into Eq. ([2|), one arrives at 
7„ = 0.67 mJ/isT^cm^ or 7o,ci ~ 41.05 mJ/K^-mol 
< 7„. Using instead Eq. ^ one obtains a slightly 
smaller value 70,01 ~ 36.95 mJ/K^-mol which however 
again is significantly smaller than our nominal value 
7„ sa 94 mJ/K^-mol. Thereby in both cases of Eqs. (|2l3p 
the experimental value Ap = 3x 10~2^J7cm has been used 
as reported in Ref.l25l . Thus, our empirical value of about 
60 mJ/K^-mol [20] can be regarded as reasonable num- 
ber. In view of the idealized electronic structure in terms 
of cylindric FSS adopted above, a more realistic and so- 
phisticated multiband analysis is desirable. To illustrate 
this point we consider the simple case when all four FSS 
would give the same contribution to the resistivity and to 
the specific heat. Applying the Kadowaki- Woods relation 
first to such a hypothetical single FSS we would arrive 



finally at 73.9 mJ/K^-mol in case of Eq. ([3]). Since differ- 
ent individual residual resistivities pii = l — A lower these 
values, we regard these two numbers as upper and lower 
bounds for a more realistic 700I somewhere in between. 
More theoretical microscopic studies including also the 
determination of the individual residual resistivities are 
necessary to improve the accuracy of these Kadowaki- 
Woods type relations for pronounced multband systems. 
Note that our empirical value of about 60 is very close to 
the mean value 0.5(36.95-1-73.9)= 55.45 and 78.57 in case 
of Eq. ©. A more detailed consideration will be given 
elsewhere. Anyhow, considering also available data for 
the in-plane penetration depth (or the condensate den- 
sity) in the superconducting state at very low T, one 
arrives at similar estimates (sec below) and wc strongly 
believe that the nominal value of about 94 mJ/K^/mol 
given above and similar numbers found in the recent lit- 
erature [2l[ as well, do significantly overestimate the con- 
tribution from the itinerant electrons which bear the su- 
perconductivity. 

The nominal value 7„ should be compared with the 
calculated quantity of 7b= 10.2 to 13.0 mJ/mol from 
DFT (densit y f unctional theory) based band structure 
calculations [251 ? , regarded as the unrenormalized 
bare quantity. The renormalization happens in two steps 
at different energy scales: a first one at high-energies is 
governed by the Coulomb interaction and/or Hund's rule 
coupling and a second one at low-energies governed by 
the interaction of the quasi-particles with various bosonic 
excitations (phonons, paramagnons, magnons etc.). The 
high-energy renormalization yields for typical transition 
metals a mass enhancement by a factor two to three as 
evidenced by a general band squeezing as observed for 
instance in ARPES measurements 0] or in optical 
measurements comparing calculated and measured un- 
screened plasma frequencies. In fact, taking a typical 122 
experimental in-plane plasma frequency of 1.55 eV [28| to 
be compared with the calculated DFT value of 2.56 eV to 
2.58 eV (25], [2^ mentioned above. Then, this high-energy 
mass enhancement can be estimated by a large factor of 
77 « 2.7 in accord with the high-energy band "squeezing" 
factor of about 2 to 3 as seen by ARPES 

Thus, one 

is left with an effective quasi-particle (qp) 7qp quantity 
of about 30 to 40 mJ/mol to be compared with our 
emprical estimate of about 60 mJ/mol K^: 



7ol = 7qp[l + Aph + Asf], 
7qp = Vlb , 

V ~ ^pl,DFT/^pl,opt ' 



(4) 



where Aph is the electron- phonon coupling constant, Agf is 
the enhancement due to spin fluctuations (antiferomag- 
netic paramagnons), and 77 > 1 denotes the high-energy 
renormalization. In case of Fe-based superconductors the 
conventional electron-phonon interaction is weak yielding 
Aph < 0.2 |3l| which is insufficient to explain the large 701- 
value obtained from specific heat measurements. A sim- 
ilar value has been found also for KFe2As2: Aph ~ 0.17 
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(details of this DFT based calculation will be given else- 
where). Taking this into account, we may finally estimate 

that Asf ^ 1 in the case of KFe2As2. 



2. Penetration depth and condensate density 

The conclusion about weak electron-boson coupling is 
also supported using the experimental value of the in- 
plane penetration depth extrapolated to T = 0: Aab,L ~ 
203 nm (measured at 50 mK) [s^]- Following Refs. 
[27l . Issl . [34 one has from the renormalized plasma fre- 
quency which enters the penetration depth rewritten in 
convenient units 



^^pi[cV]AL[nm] = 197. DNZ,,, > 1 , (5) 

where N = ntot is the reciprocal number of the con- 
duction electron density involved in the superconducting 
condensate, Z^^ « (1 -I- Atot(O)) describes the dynami- 
cal mass renormalization and D — (1-1- 5) (14-/) with 
S, f > 1 describes the effect of disorder and fluctuations 
of competing phases. In the clean limit one has 6 0. 
For the sake of simplicity we will ignore the influence 
of fluctuations. Using hil.pi ~ 1-55 eV for the expected 
unscreened experimental plasma frequency, (i.e. the ex- 
perimental high-T plasma energy with no or small renor- 
malizations due to the electron-boson couplings) one has 



2.64^ = (1 + Aph + A,f ) (1 + S) 

"-tot 



(6) 



where 5^1/3 measures the disorder and the recipro- 
cal gap amplitude related parameter close to that in the 
clean limit (i.e. 5 <C 1). Furthermore for the sake of sim- 
plicity we will assume that all electrons are involved in 
the superconducting condensate, i.e. iV = 1 (just for illus- 
tration see also the special case Us/n sa 0.74 mentioned 
in our remark (35j). Then, one arrives at the constraint 
M: 



Atot = \h + Asf w 0.97, or A^f « 0.8 



(7) 



in accord with a close estimate from 7e; ~ 60 mJ /K^ mol 
and 7h sa 10.5 mJ/K^mol from DFT calculations for the 
bare density of states (26| . 



D. The upper critical fields Hc2(T) and their 
anisotropy 

Figs. 5-7 summarize the T-dependence of the specific 
heat data Cp of the investigated KFe2As2 single crystal for 
various magnetic fields applied parallel and perpendicu- 
lar to the a6-planc. With increasing applied magnetic 
field in both directions, the superconducting anomaly 
shifts and broadens systematically to lower T and is also 
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FIG. 5: (Color) The electronic specific heat coefficient Cei/T 
of KFe2As2 (after subtracting the phonon contribution) for 
both directions -B || c and B \\ ab a& shown in (a) and (b), 
respectively. In order to determine the Tc of KFe2As2, an 
entropy - conserving construction has been used as shown 
with a green line in the inset of (b). The insets of the upper 
and lower panel show two data sets with the same Tc value 
for the two directions confirming our anisotropy ratio F ~ 5. 



reduced in height. In an applied magnetic field of 9 T, su- 
perconductivity is completely suppressed for both direc- 
tions of our crystal. In order to to analyze the phase di- 
agram of the field dependence of Tc, we used an entropy- 
conserving construction of the electronic specific heat to 
determine Tc of both orientations as shown in Fig. 5. 
Then, in a very first step the upper critical field and its 
slope near Tc can be estimated by the Ginzburg Landau 
(GL) equation [13, HI] (strictly speaking valid near T^, 
only): 



Hc2(0) 



(8) 



where t — T/Tc- The upper critical field values at T = 



have been evaluated to HqH'^ (0) 



1.8 T and Mo^cf^O) 



— 8.6 T and the fits are shown via dashed black lines in 
Fig. 6. In the case of i7c2 obtained from ac susceptibil- 
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ity data (see above Fig. [T]) we derive at slightly different 

values of Ato-ffc2^(0) w 1.9 T and ^loH[f\Q) w 7.5 T, re- 
spectively (dashed red curves). 

It is interesting to compare the obtained and extrap- 
olated to T = anisotropy ratio for the upper critical 
field of 4 to 5 with that ratio for the penetration depth 
[39| (i.e. \L,ab = 194, 3 nm and \l,c =510.3 nm taken at 
T =20 mK) which yields 2.63, only. In fact, for a simple 
one-band or separable multiband models [12] including a 
phenomenological mass anistropy, one would expect 



mab 



1/2 



Al.c(O) ^ gc2||afc(0) 
AL,afc(0) H,2\\M 



(9) 



From full relativistic DFT-calculations an out of plane 
plasma frequency of 0.61 eV has been obtained [2^ which 
suggests an mass anisotropy of 4.38 slightly exceeding 
the value of 3.27 for Ba-122 [13 Thus, the observed 




T(K) 



FIG. 6: (Color online) Phase diagram of fioHc2{T) for 
KFe2As2 with the magnetic field apphed parallel and perpen- 
dicular to the c axis from specific heat (black symbols) and ac 
susceptibility measurements (red symbols). Solid lines: the- 
oretical curves based on the WHH model (a = and Xso ~ 
0.1). Dashed curves: GL theory. Inset: T-dependence of _ffirr 
obtained from ac susceptibility measurements as discussed in 
the text. 



anisotropy derived from the upper critical fields exceeds 
this value whereas the penetration depth gives a slightly 
smaller value. We ascribe this small deviation of our em- 
pirical F from the simple mass anisotropy to 

(i) the anisotropy of the pairing interaction and conse- 
quently also of the order parameter and /or oppositely 
of the depairing interaction, in Ref. |39| which might 
additionally enhance Hc2\\ab and suppress XL,ab or vice 
versa the corresponding c-components. For instances, the 
anisotropic screening and significantly anisotropic plasma 
frequencies might cause an anisotropic Coulomb pseu- 
dopotential /x* . The in-plane anisotropy observed for or- 
dered magnetic structures should in case of a magnetic 
spin fluctuation bases mechanism act in a similar way. 

(ii) Furthermore, one should take into account, that 
strictly speaking, the upper critical fields and the pene- 
tration depth at r = probe various subgroups of elec- 
trons with different Fermi-velocity dependent weights: 
whereas the penetration depth probes more sensitively 
fast electrons A^^^ oc ^pi^i, where flpi^i denotes the cor- 
responding i*^ subgroup plasma frequency and the total 
penetration depth is given by A^^ = X^iA^^j. In con- 
trast, the upper critical fields are more sensitive to slow 
electrons since ffc2||ab {^o/vxVz), where $0 denotes 
the flux quantum. Finally, 

(iii) anisotropic impurity scattering rates might also 
affect Fq. 

Another possibility to estimate roughly the upper criti- 
cal field -ffc2(0) is to consider the single-band Werthamer- 
Helfand-Hohenberg (WHH) formula [4l| with the Maki 
parameter a = and Xso = 0.1 which is a reasonable 
value of spin orbital scattering for Fe-based superconduc- 
tors [l^l. As shown with solid lines in Fig. [HI the specific 
heat and ac magnetization Hc2 data for B \\ ab arc per- 
fectly described by the WHH model with an average slope 
-d{fj,oH^^''^)/dT « 2.8(2) T/K while for H \\ c the single- 



band WHH model with 



-d{fioH^J)/dT = 0.55(5) T/K 
underestimates the experimental data obtained from the 
specific heat measurements (see lower panel of Fig. (5]). 
This deviation is particularly visible in the case of Hc2 
obtained from the ac magnetization data where we used 
-d{fiQH^^f)/dT =0.6(1) T/K. From these values the up- 
per critical fields /Xo^^c2(0) are found to be ~ 1.4 T and 
~ 7 T for the c and ab direction, respectively. The ob- 
served small difference between Hc2 obtained from the 
specific heat and the ac magnetization data is expected 
since these methods naturally implies different criteria 
for Tc determination. 

In general, in case of multi-band superconductivity 
the low-T i/f^-curve may exceed the single-band WHH 
predictions [43]. Therefore, we suppose that the ob- 
served deviation from the single band WHH model is 
related to multi-band effects. Additionally, indication 
for a two-band like behavior of our single crystal was ob- 
served in zero field specific heat measurements (see be- 
low). Anyhow, using typical renormalized Fermi veloci- 
ties vp ^ A X 10*m/s derived from preliminary ARPES- 
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data [6[ and Tc = 3.5 K, one estimates also, in princi- 
ple, within a two-band approach adopting s-symmetry 
[il . [45} , a slope- value 



E. Aspects of the electronic specific heat in the 
superconducting state: the residual linear specific 
heat and the jump at Tc 



(10) 



where c\ — )■ C2 — )■ 1/2 and vp ^ ■\/(2)ui, y/j2)v2 in the 
case of a dominant interband pairing and C(3) sa 1.202, 
resulting in -dH^2/dT=Qm T/K near Tc which is al- 
ready very close to our experimentally determined value 
and is also in accord with the renormalized Fermi velocity 
of 4.-10^cm/s using the total bare velocity 1.77-lO^cm/s 
from the full relativistic (not spin polarized) LDA cal- 
culations and the FSS averaged renormalizations con- 
tained in the intrinsic 7orvalue of about 60 mJ/K^mol 
estimated above. In comparison, the reported values de- 
termined via detailed resistivity studies on KFe2As2 sin- 
gle crystals yield lower values, i.e., = 1-25 T and 
H'g = 4.47 T, where a low value of = 2.8 K, has been 
reported Q. The anisotropy of the slopes near Tc as 
measured of about 5.35 is very close to the value found 
here: 5.09. The reported larger value of 6.8 seems to 
be a consequence of the extremely high anisotropic spin- 
orbit coupling = 0.36 for B \\ ah and 00 for B \\ c 
adopted in Ref. |3 in analyzing their dataji^. The re- 
ported larger absolute slope values might be interpreted 
as a hint for an impurity driven transition to an s-wave 
superconductor with (A)fs 7^ with pair-breaking (see 
Eq. (A3) in Ref. [T2h . From our studies, further informa- 
tion about the anisotropy of KFe2As2 single crystals can 
be obtained, which is F = i?c2/^c2 ^ (see also the 
insets of Fig. 5). Surprisingly, this anisotropy value is 
comparable with F- values of e.g. NdFeAsO0.82F0.i8 [13 
and LaFcPO 0| showing a more anisotropic electronic 
structure (LaFeAsO: 9.2 to 10.8 and LaFcPO: 4.16 to 5.04 
and might be therefore ascribed to opposite anisotropics 
of the order parameter. On the other hand, it is consid- 
erably larger than a typical value of F ^ 2 and 2.6 found 
for nearly optimally hole doped BaKFe2As2 [igl.lsoj. but 
lower than the ones determined for SmFeAsO0.85F0.15 
and La(0,F)FeAs thin films ^EM- 

The T-dependence of the irreversibility field Hirr ob- 
tained from x"v are shown in the inset of Fig. 6 (see 
above). The low value of i/f/^ for _B || ab is related 
with a large anisotropy and a weak pining as expected in 
the case of clean single crystals. We attribute the Hlj^j. 
with a peak effect in the T-dependence of the critical cur- 
rent Jc for iJ II c in accord with similar observations on 
YBCO single crystals [l^ which exhibit a rather similar 
anisotropy of the upper critical field and therefore also a 
similar pinning behavior can be expected. 



The jump height of Acoi/Tc « 45.6 mJ/mol at 
Tc is found from our zero-field electronic specific heat 
data. This value exceeds the value which has been re- 
ported for a polycrystalline KFe2As2 sample [H but by 
a factor of 2 lower than the one obtained for the nearly 
optimally- hole doped Bao.6Ko.4Fe2As2 [i^. For our es- 
timated 7oi ^ 60 mJ/mol-K^ the ratio Acei/jeiTc found 
to be enhanced as compared with the use of the nominal 
value: near about 0.76 vs. 0.49 (see Fig. 4), but it is still 
significantly lower than the result of the BCS weak cou- 
pling approximation; Acei/7eiTc=1.43 [s^ . which points 
towards a multiple band (gap) with s,p- or c?-wave na- 
ture of superconductivity. In particular, it is close to the 
value reported for the p-wave superconductor SrRu04 
which exhibits 0.73 [ssj . 

In a clean situation with negligible pair breaking cf- 




FIG. 7: (Color online) The T-dependence of the specific heat 
of KFe2 As2 in various applied magnetic fields up to 9 T par- 
allel to the c axis (a) and parallel to the ab plane (b). The 
insets of the upper and lower panel show a zoom into the 
superconducting state for both directions. 
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fects, the reduced jump of the specific heat Acd /Tc 701 
as compared with that of a single s-wave superconduc- 
tor might be related to unconventional superconductiv- 
ity with nodes as discussed above and/or to pronounced 
muhiband character with rather different partial densi- 
ties of states and gaps. Furthermore, in relative dirty 
systems the unconventional superconductivity might be 
driven into an s-wave state. To illustrate the multi- 
band character we adopt here for the sake of simplic- 
ity a simple effective weak coupling s-wave model like in 
Ref. [23. Another interesting issue we would like to ad- 
dress is what happened with the "extrinsic" linear spe- 
cific heat at very low T . Thus, fitting the electronic part 
of the specific heat within a two-band model (see the 
blue curves in Fig. 8) admitting also a "residual" linear 
Sommerfeld part, we arrive at a relative large value of 
lrcs{T — > 0) ~ 15 mJ/K^ mol which might be related to 
(i) an "extrinsic" pair breaking contribution somewhat 
suppressed deep in the superconducting state [35l |. Any- 
how, we admit that the adopted s-wave analysis might 
provide only an upper limit, since for an unconventional 
pairing symmetry the spectral weight at low-T is en- 
hanced. Furthermore the final density of states intro- 
duced by pair-breaking induced subgap states might also 
contribute to such a value. Thus, more sophisticated 
multiband models including interacting pair-breaking im- 
purity states are necessary to settle this highly interest- 
ing problem. Due to its complexity it is however far be- 
yond the present paper. Specific heat measurement below 
0.2 K might be helpful to refine the value of ^res- In this 
contect the observation of substantial residual terms in 
other pnictide or chalcogenide superconductors is note- 
worthy. For instance in the systems FeTcg. 57809.45 and 
Co-doped Ba-122 also a relatively large (8% and 25% , 
respectivel y) r esidual linear contribution have been ob- 
served [53.l55j. 

Finally, for completeness we list the gap values ob- 
tained in the present simple model for analyzing the T- 
dependence of our zero-field specific heat measurements 
down to 400 mK. The normalized zero-field electronic 
specific heat Cd/jnT is shown in Fig. 8. First, we com- 
pare our data to the single-gap BCS theory (i.e. the weak 
coupling approach using Ao/ksTc = 1.76 at Tc) and find 
that a single BCS-gap cannot be reconciled with our ex- 
perimental data. From the plot in Fig. 8 it is evident that 
the specific heat jump is much lower than the BCS weak- 
coupling value for an s-wave superconductor. Below Tc, 
systematic deviations of the theoretical curve from the 
experimental data have been observed at both low-T 
and around Tc- This clearly indicates that Cei/jeiT of 
KFe2As2 cannot be described by the weak coupling BCS 
(single-)gap. 

Since a single-gap scenario cannot describe our data, 
we applied a phenomenological two-gap model in line 
with multi-gap superconductivity in many compounds of 
the FeAs family reported by various experimental and 
theoretical studies on different compounds within this 
family (2|. [56l - [58| . We have analyzed our data utilizing the 



generalized a-modcl which has been proposed to account 
for the thermodynamic properties in multi-band, multi- 
gap superconductors like e.g. MgB2 We remind the 
reader that in this approach the one-band expression: 

-^--^r^ / [/ln/ + (l-/)ln(l-/)]dy, (U) 



is straightforwardly generalized to the two-band case and 
entropy conservation is adopted for each band. In Eq. 
([TT|) the Fermi- function is denoted by / = [exp(/3£^ -|- 
1)]^^, /3 = (kBT)^^ and the energy of the quasiparti- 
cles is given hy E = [e^ + A'^{t)]^-^ with e being the 
energy of the normal electrons relative to the Fermi sur- 
face. The integration variable is y = c/Aq. (S) and 
(C) is the thermodynamic properties and t = T/Tc is 
the reduced temperature. In Eq. (fTTl the scaled gap 
a ~ Ao/fcsT is the only adjustable fitting parameter. 
The temperature dependence of the gap is determined 
by A(t) = AoS{t), where S{t) is approximately described 
by the data taken from the table in Ref. In case of 
two gaps the thermodynamic properties are obtained as 
the sum of the contributions from the two gaps, i.e., ai = 
Ai{0) /ksTc and a2 = A2{0) /ksTc with their respective 
weights 71 /7m and 72/70I. 

To calculate the theoretical curves Cd/jciT the param- 
eters Ai, A2, their respective ratios 71 and 72 and the 
ratio "fres/jei are left for free in fitting {jres represents 
the non-negligible residual value at low-T). The best 
description of the experimental data is obtained using 
values of Ai/fesT^ = 0.46 and Az/fcsTc = 1.75. The 
calculated specific heat data are represented by the solid 
blue line in Fig. 8 (upper panel. Anyhow, small rela- 
tive jumps are not compatible with the strong coupling 
scenario estimated in Fig. 4 for the case of no extrinsic 
contributions. Therefore we performed a second analysis 
where the effective extrinsic linear contibution necessary 
for a weak coupling scenario has been subtracted fromn 
the row data. The result is shown in the lower panel of 
Fig. 8. Then both gaps do slightly increase: 1.8 K and 
6.2 K. 

Naturally, the obtained gap values are smaller than 
the largest gap in e.g. the optimally hole doped 
Ban.fiKn.4Fe2As2 as observed by ARPES investiga- 
tions [6l| but comparable to the two-band s-wave fit for 
the penetration depth data (H |j c: 1.28 K and 5.57 K 
but do not clearly exceed the corresponding values for the 
isomorphic compound RbFe2As2 system with a slig htly 
lower Tc-value of 2.52 K, only: 1.74 K and 5.7 K. 
In our opinion this might refiect the presence of nodes in 
the superconducting order parameter of KFe2As2. Any- 
how, a detailed comparison of these two closely related 
systems would be very interesting, especially, if in fact 
it would be confirmed that the symmetry of the order 
parameter would be different. 
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FIG. 8: (Color online) Upper panel: Fit under the assump- 
tion of no extrinsic contribution to the linear specific heat 
(see also Figs. 3 and 4). The normalized superconducting 
electronic specific heat Cci/i^eiT) of KFe2As2 as a function 
of the reduced temperature t = T /Tc- The red line repre- 
sents the theoretical curve for the single-band weak coupling 
BCS case (Ao/fcsTc = 1.76). The blue line shows the curve 
of the nodeless weakly coupled two-gap model fit; for details 
see text. The inset shows our electronic specific heat data in 
comparison with data by Fukazawa et at, [2l|] together with 
an another sample which has been shown that our Tc is sim- 
ilar to the investigated sample. Lower panel: fit assuming 
that a significant contribution to the linear specific heat is 
not intrinsic: using e.g. 7ci = 60 mj /mol taken from Ref. 
HO (see Figs. 3 and 4.). 

Although a clear picture is still missing for the case of 
KFe2As2, it is important to emphasize that our system 
definitely underlies multiband superconductivity, prob- 
ably in the weak coupling regime. However, from spe- 
cific heat data alone it is difficult to be sure whether 
nodes exist or not, since in the case of multi-band super- 
conductivity low-energy quasi-particle excitations can be 
always explained by the contribution from an electron 
group with a small gap. We believe that further experi- 
mental studies such as specific heat well below 400 mK, 
ARPES and transport investigations at very low T, will 
be helpful to elucidate the nature of superconductivity in 



KFesAss. 

IV. CONCLUSIONS 

In summary, KFe2As2 was investigated by ac suscep- 
tibility and low-T specific heat measurements on high- 
quality single crystals grown by a self-flux technique. 
The specific heat jump was found to be Acd/Tc ^ 
45.9 mJ/mol and the nominal Sommerfeld coefficient 
7n = 94(3) mJ/mol K^. However, several theoretical 
considerations including two recently proposed modified 
Kadowaki- Woods relations as well as the observation of 
a significant linear in T residual term point to a sig- 
nificantly smaller value for the itinerant quasi-particles 
of about 60 mJ/K^ mol. This way the strongly corre- 
lated " heavy- fermion-like" scenario suggested for K-122 
in the literature should be revisited. In this context the 
elucidation of the "external" subsystem responsible for 
that difference is a challenging problem to be considered 
elsewhere. The total electron-boson coupling constant 
Atot = Aph + Asf ^ 1 avaraged over all Fermi surfaces 
excludes strong coupling. The calculated weak electron- 
phonon coupling of about 0.17 points to a dominant spin- 
fluctuation mechanism and unconventional superconduc- 
tivity. 

The magnetic phase diagram has been studied yield- 
ing values for the upper critical fields /ioi?c2(0) ~ 1-4T 
and /io-f^c2'(0) « 7T for the c axis and ab plane, respec- 
tively. The resulting anisotropy of KFe2As2 near Tc lies 
around F = H^^'^ / H^!^ ^ 5 which slightly exceeds the 
mass anisotropy as derived from DFT-clectronic struc- 
ture calculations as well as the anisotropy of the pen- 
etration depth. Additionally, the T-dependence of our 
zero-field electronic specific heat c^i cannot be described 
within single-band weak-coupling BCS theory. 

For a full understanding of the gap structure of 
KFe2As2 as well as of the high values of Ce\/jeiT at low 
temperatures, further specific heat measurements at very 
low T < 400 mK and/or low-T ARPES and transport 
studies will be helpful. Finally, the irreversibility field 
Hirr derived from ac susceptibility data has been investi- 
gated. The double-maximum in x"viT) for H\\c suggests 
the presence of a peak effect in the T-dependence of the 
critical current. 
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